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TITLE OF THE INVENTION 
ANALYTICAL MODEL CONVERSION METHOD 

FIELD OF THE INVENTION 
5 The present invention relates to a method of 

creating a shell element analytical model for 
two-dimensional analysis from a three-dimensional solid 
analytical model by a finite element method on the 
basis of three-dimensional geometric data designed by 
10 three-dimensional CAD. 

BACKGROUND OF THE INVENTION 

As three-dimensional CAD is recently becoming 
popular, preparation of a quick analytical model 

15 creating method using three-dimensional CAD geometric 
data has been demanded. As is represented by linear 
structure analysis, a three-dimensional solid 
analytical model can be created by tetrahedral element 
division by a finite element method based on 

20 three-dimensional CAD geometric data so that analysis 
can easily be executed. 

However, if thin-walled structures or exterior 
components are to be analyzed by using 

three-dimensional solid analytical models, as in, e.g., 
25 plastic injection molding CAE 

(filling/packing/cooling/warp analysis program), the 
temperature distribution in the direction of plate 
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thickness must be accurately calculated. To do this^ 
element division in the direction of plate thickness in 
about five to 10 layers is necessary. At this time, 
the total number of elements of analytical models is 
5 more than 5,000,000. It is not practical in terms of 
calculation cost from the viewpoint of current 
computing capability. In many cases, therefore, a 
neutral plane is created manually from a 
three-dimensional solid model. Calculation is 
10 performed by using an analytical model obtained by 

dividing the neutral plane by triangular or rectangular 
two-dimensional shell elements. 

Examples of prior arts of the neutral plane and 
neutral plane element generation method are as follows. 
15 1) Method of generating a neutral plane in the 

direction of plane thickness of a three-dimensional 
solid geometric model 

Two-dimensional neutral plane surface data is 
manually generated for two planes that constitute the 
20 plane thickness of geometric data from 

three-dimensional solid CAD data. A neutral plane can 
be generated by using a CAD/CAE tool represented by, 
e.g., IDEAS, PATRAN, and FEMAP. There also exists a 
method of generating a two-dimensional shell element 
25 model for the neutral plane shape. 

2) Method of generating a hexahedral solid 
element and generating a two-dimensional shell element 
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by using an intermediate node 

In a method of prior art, a hexahedral solid 
element is generated from three-dimensional solid CAD 
data by using a preprocessor such as a finite element 
5 method, and a shell element is generated by using an 
intermediate node. 

However, the above -described neutral plane and 
neutral plane element generation methods have the 
following problems. 
10 1) Method of generating two-dimensional neutral 

plane surface data from three-dimensional solid CAD 
data 

It is difficult to accurately generate a neutral 
plane for a portion where planes cross each other, a 
15 portion where the plate thickness changes, or a 

thick-walled portion. In addition, no plate thickness 
information can be added to the generated neutral plane 
surface even when the plate thickness changes in a 
tapered manner. Hence, plate thickness information 
20 must be manually input after element division. 

2) Method of generating a hexahedral solid 
element and generating a two-dimensional shell element 
by using an intermediate node 

A hexahedral solid element is hard to 
25 automatically generate even by using a preprocessor. A 
hexahedral solid element needs to be manually 
generated. Model creation is therefore very 
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t ime - consuming . 

As described above, in both the method of 
generating a neutral plane in the direction of plane 
thickness of a three-dimensional solid geometric model 
5 and generating a two-dimensional shell element model 
for the neutral plane shape and the method of 
generating a hexahedral solid element in a 
three-dimensional solid geometric model and generating 
a two-dimensional shell element by using an 
10 intermediate node, analytical model creation is very 
time-consuming . 

SUMMARY OF THE INVENTION 
The present invention relates to a method of 
15 automatically creating a two-dimensional shell element 
analytical model on the basis of three-dimensional CAD 
geometric data. It is therefore a feature of the 
present invention to largely shorten the time for 
creation of a two-dimensional shell element analytical 
20 model from three-dimensional CAD geometric data. 

According to the present invention, there is 
provided an analytical model conversion method of 
converting a three-dimensional analytical model into a 
two-dimensional analytical model, comprising generating 
25 tetrahedral solid elements for an input 

three-dimensional geometric model, and connecting 
intermediate nodes of sides that extend in a direction 
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Of plate thickness in each tetrahedral solid element to 
generate a triangular or rectangular shell element. 

Other features and advantages of the present 
invention will be apparent from the following 
5 description taken in conjunction with the accompanying 
drawings, in which like reference characters designate 
the same or similar parts throughout the figures 
thereof . 



10 BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a flow chart showing the outline of a 
neutral plane element generation method; 

Fig. 2 is a view showing a case wherein a 
triangular neutral plane element is generated in a 
15 tetrahedral element; 

Fig. 3 is a view showing a case wherein a 
rectangular neutral plane element is generated in a 
tetrahedral element ; 

Fig. 4 is a view showing an example of 
20 tetrahedral element division of a shape having a plate 
thickness difference; 

Fig. 5 is a view showing another example of 
tetrahedral element division of a shape having a plate 
thickness difference; 
25 Fig. 6 is a view showing an example of 

tetrahedral element division of a rib shape; 

Fig. 7 is a view showing another example of 
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tetrahedral element division of a rib shape; 

Fig. 8 is a view showing still another example of 
tetrahedral element division of a rib shape; 

Fig. 9 is a view showing an example of 
5 tetrahedral element division of an L shape; 

Fig. 10 is a view showing another excunple of 
tetrahedral element division of an L shape; 

Fig. 11 is a view showing an example of 
tetraJtiedral element division of a cross shape; 
10 Fig. 12 is a view showing an exsunple of 

tetrahedral element division of a shape having crossing 
radial portions; 

Fig. 13 is a flow chart showing details of the 
neutral plane element generation method; 
15 Fig. 14 is a view showing the three-dimensional 

CAD geometry of a component of a fixing toner container 
used in a laser beam printer (LBP)^ which is designed 
by IDEAS; 

Fig. 15 is a view showing the three-dimensional 
20 CAD geometry of a component of a fixing toner container 
used in a laser beam printer (LBP), which is designed 
by IDEAS; 

Fig. 16 is a view showing divided elements 
obtained by automatically dividing the 
25 three-dimensional CAD geometric data into tetrahedral 
elements by using a preprocessor function for finite 
element analysis of IDEAS so as to form a 
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Single -layered structure in the direction of plate 
thicknes s ; and 

Fig. 17 is a view showing a shell element model 
for which neutral plane shell elements are generated 
from the tetrahedral elements on the basis of a neutral 
plane element division method. 

Fig. 18 is a block diagram showing the 
construction of a teiminal which runs a system for 
simulating the conveyance of the medium according to 
the first embodiment. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
As described above, the present invention relates 
to a method of automatically creating a two-dimensional 
shell element analytical model on the basis of 
three-dimensional CAD geometric data. With this 
method, time for creation of a two-dimensional shell 
element analytical model from three-dimensional CAD 
geometric data can largely be shortened. The 
embodiment of the present invention will be described 
below in detail with reference to the accompanying 
drawings - 
[ Embodiment ] 

Fig. 1 shows one schematic process of an 
apparatus of the present invention. Fig. 18 is a block 
diagram showing the construction of the apparatus of an 
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embodiment of the present invention • 

A central processing unit (CPU) 101 performs the 
overall control of the terminal on the basis of 
programs expanded in a main memory 103. An input 
5 device 102 is a pointing device including a key board, 
a mouse, etc. The main memory 103 is constructed of a 
random access memory (RAM) or the like and serves as a 
work memory for, for example, expanding the programs. 
A display 104 is constructed of a cathode-ray tube 
10 (CRT) monitor, a liquid crystal display, or the like. 
An auxiliary memory 105 is constructed of a hard disk 
drive or the like and stores various programs for 
operating a server (or the terminal) and various 
databases. A communication device 106 is an interface 
15 for providing connection to a network. 

First, as shown in step SI in Fig. 1, assume that 
an object shape is designed by three-dimensional CAD so 
that three-dimensional CAD geometric data is available. 
Processing starts from a state wherein the 
20 three-dimensional geometric data is loaded by a 

preprocessor such as IDEAS available from EDS PLM 
solutions or PATRAN available from MSC, which is used 
for finite element analysis and is capable of 
tetrahedral element division. 
25 In step S2 in Fig. 1, tetrahedral elements are 

generated from the shape by using the automatic 
tetrahedral element division function of the 
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preprocessor. At this time, the length of one side of 
each tetrahedral element must be designated such that 
the shape is divided in the direction of plate 
thickness to form a single -layered structure. This can 
5 be done by designating the maximum plate thickness of 
the shape. The node coordinate information and element 
constituent node information of the tetrahedral 
elements generated by element division are written in 
external files. If the shape has a thick-walled 
10 portion, and the number of layers in the direction of 
plate thickness should be two or more, the shape is 
corrected in advance on the three-dimensional CAD side. 

In step S3 in Fig. 1, on the basis of the node 
coordinate information and element constituent node 
15 information of the tetrahedral elements generated in 
step S2, an intermediate node is generated in each of 
sides that constitute a tetrahedral element. This 
processing may be executed even by designating 
secondary tetrahedral elements (tetrahedrons having an 
20 intermediate node in each side) in automatic 
tetrahedral element division in step S2. 

In step S4 in Fig. 1. triangular or rectangular 
neutral plane shell elements are generated in the 
direction of plate thickness of the shape by using the 
25 intermediate node information in step S3. VRien the 

tetrahedral elements are generated such that the number 
of layers formed by dividing the shape in the direction 
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Of plate thickness becomes one, a tetrahedral element 
generated between the upper and lower planes in the 
direction of plate thickness can basically take two 
forms. As the first form, a plane and a corresponding 
apex of the tetrahedral element are located on the 
upper and lower planes of the shape, as shown in 
Fig. 2. As the second form, two sides of the 
tetrahedral element are located on the upper and lower 
planes of the shape, as shown in Fig. 3. The neutral 
plane shell element is triangular in the former case 
and rectangular in the latter case. 

In step S5 in Fig. 1, for each of the triangular 
or rectangular neutral plane shell elements generated 
in step S4, the plate thickness until the geometric 
surface in the direction of normal of the neutral plane 
shell element is calculated and defined as the plate 
thickness of the shell element. 

Finally, in step S6 in Fig. 1. analysis input 
data is created by adding boundary conditions or 
analysis conditions corresponding to the type of 
analysis to the node coordinates, element constituent 
node, and element plate thickness information of the 
triangular or rectangular shell elements generated in 
accordance with the above -described procedures. Then, 
analysis is executed. 

The basic outline until analysis input data is 
created by generating neutral plane shell elements from 
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a three-dimensional shape has been described above. 

The processing must cope with the geometric 
elements of an actual complex product shape. A method 
for it will be described in detail. 

Tetrahedral element generation patterns and 
neutral plane shell element generation methods will be 
described in consideration of geometric elements having 
plate thickness differences or rib structures of 
various types. 

Figs. 4 and 5 show tetrahedral element generation 
patterns when a shape has a plate thickness difference 
on the geometric section in the direction of plate 
thickness. Basically, each corner portion of the shape 
has a tetrahedral element whose two or more element 
planes are located on the outer surface of the shape, 
as shown in Figs. 4 and 5. This element will be 
referred to as a "corner element". The "corner 
element" (Fig. 5) at the plate step portion can be 
excluded. However, a "corner element" at an end 
portion of the shape cannot be omitted. Hence, the two 
types of corner elements must be distinguished. 

In addition, a tetrahedral element (Fig. 5) whose 
element planes are not located on the outer surface of 
the shape at all is formed. This element will be 
referred to as an "internal element". To generate a 
neutral plane shell element in this "internal element", 
the constituent nodes of the neutral plane shell 
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elements in adjacent tetrahedral elements and their 
continuity must be taken into consideration. 

Figs. 6. 7, and 8 show tetreihedral element 
generation patterns when a shape has a T-shaped section 
5 because of addition of a rib. Even in these examples, 
a tetrahedral element as an "internal element" whose 
element planes are not located on the outer surface of 
the shape at all is generated. Hence, a shell element 
must be generated in consideration of the constituent 
10 nodes of the neutral plane shell elements in adjacent 
tetrahedral elements and their continuity. 

Figs. 9 and 10 show tetrahedral element 
generation patterns when a shape has an L- shaped 
section. Even in these examples, "corner elements" and 
15 "internal elements" are generated. 

Figs. 11 and 12 show tetrahedral element 
generation patterns when a shape has a cross shape or 
crossing radial portions. For these examples, 
basically, that a plurality of "internal elements" are 
20 generated must be taken into consideration. 

As described above, to generate neutral plane 
shell elements by generating tetrahedral elements in a 
three-dimensional shape such that the number of divided 
layers in the direction of plate thickness becomes one, 
25 the neutral plane shell elements are generated by 

making remark on the above -described "corner elements" 
and "internal elements" . This method will be described 
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next In detail. 

Fig. 13 is a flow chart for explaining the 
above -described flow shown in Fig. 1 in more detail. 
Fig. 13 shows a detailed flow to generate a neutral 
plane shell element from a generated tetrahedral 
element . 

First, as shown in step SIO in Fig. 13. assume 
that an object shape is designed by three-dimensional 
CAD so that three-dimensional CAD geometric data is 
available. As described above, processing starts from 
a state wherein the three-dimensional geometric data is 
loaded by a preprocessor such as IDEAS available from 
EDS PLM solutions or PATRAN available from MSC, which 
is used for finite element analysis and is capable of 
tetrahedral element division. 

In step Sll, tetrahedral elements are generated 
from the shape by using the automatic tetrahedral 
element division function of the preprocessor. At this 
time, the length of one side of each tetrahedral 
element must be designated such that the shape is 
divided in the direction of plate thickness to form a 
single- layered structure. This can be done by 
designating the maximum plate thickness of the shape. 
The node coordinate information and element constituent 
node information of the tetrahedral elements generated 
by element division are written in external files. If 
the shape has a thick-walled portion, and the number of 
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layers In the direction of plate thickness should be 
two or more, the shape is corrected in advance on the 
three-dimensional CAD side. 

In step S12, on the basis of the node coordinate 
information and element constituent node information of 
the tetrahedral elements generated in step Sll, an 
intermediate node is generated in each of sides that 
constitute a tetrahedral element. This processing may 
be executed even by designating secondary tetrahedral 
elements (tetrahedrons having an intermediate node in 
each side) in automatic tetrahedral element division in 
step Sll. 

In step S13. for all the divided tetrahedral 
elements, a table of the numbers of elements that are 
adjacent to the planes of the tetrahedral elements is 
created. In step S14. for all the divided tetrahedral 
elements, a table of planes that are not adjacent to 
the planes of the tetrahedral elements (i.e.. planes 
that constitute the surfaces of the shape) is created. 
These tables can easily be created. 

In step S15, a "corner element" as a tetrahedral 
element whose two or more element planes are located on 
the outer surface of the shape is detected by 
determining each element of the shape. In step S16, an 
"internal element" as a tetrahedral element whose 
element planes are not located on the outer surface of 
the shape at all is detected. The "corner elements" 
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and "internal elements" can also easily be detected. 

In step S17, "corner elements" that are 
unnecessary for neutral plane shell element generation 
in the geometric data divided into the tetrahedral 
elements are excluded. As already described above in 
the example of the plate step portion shown in Fig. 5, 
the "corner element" at the plate step portion is 
excluded. However, a "corner element" at an end 
portion of the shape cannot be omitted. Hence, the two 
types of corner elements must be distinguished. To do 
this, a "corner element" adjacent to an "internal 
element" is extracted first. An element to be excluded 
can be decided by determining whether two or more 
element planes of the extracted "corner element" belong 
to the geometric surface constituent element table 
created in step S14. In step S18, the element is 
excluded, and the adjacent element information of 
"internal elements" is updated. 

After this preparation, in step SI 9. triangular 
or rectangular neutral plane shell elements are 
generated in the direction of plate thickness of the 
shape by using the intermediate node information in 
step S12. When the tetrahedral elements are generated 
such that the number of layers formed by dividing the 
shape in the direction of plate thickness becomes one. 
a tetrahedral element generated between the upper and 
lower planes in the direction of plate thickness can 
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basically take two forms. As the first form, a plane 
and a corresponding apex of the tetrahedral element are 
located on the upper and lower planes of the shape, as 
shown in Fig. 2. As the second form, two sides of the 
tetrahedral element are located on the upper and lower 
planes of the shape, as shown in Fig. 3. The neutral 
plane shell element is triangular in the former case 
and rectangular in the latter case. 

In step S20, for each of the triangular or 
rectangular neutral plane shell elements generated in 
step S19, the plate thickness until the geometric 
surface in the direction of normal of the neutral plane 
shell element is calculated and defined as the plate 
thickness of the shell element. In step S21, it is 
checked whether, in all the triangles generated in step 
S19, two adjacent triangles can be converted into a 
rectangle, by considering the internal angles on both 
sides of the adjacent sides of the two triangles. If 
the two triangles can be converted into a rectangle, 
conversion is performed. 

After the above operation, analysis input data is 
created by adding boundary conditions or analysis 
conditions corresponding to the type of analysis to the 
node coordinates, element constituent node, and element 
plate thickness information of the triangular or 
rectangular shell elements generated in accordance with 
the above -described procedures. Then, analysis is 
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executed . 
[ Example 1 

An example in which the above-described 
embodiment is applied to an actual product shape will 

5 be described next. Plastic injection molding CAE 

(filling/packing/cooling/warp analysis program) will be 
exemplified here, in which a thin-walled plastic 
structure exterior component is modeled by shell 
elements and analyzed. 

10 Figs. 14 emd 15 show the three-dimensional CAD 

geometric data of a component of a fixing toner 
container used in a laser beam printer (LBP). which is 
designed by IDEAS. The basic plate thickness of this 
component is 2.5 mm. The component has a complex rib 

15 shape. 

Fig. 16 is a view showing divided elements 
obtained by automatically dividing the 
three-dimensional CAD geometric data into tetrahedral 
elements by using a preprocessor function for finite 
20 element analysis of IDEAS so as to form a 

single -layered structure in the direction of plate 
thickness. This view of divided elements is output to 
a universal file (text data) as the intermediate format 
file of IDEAS. 

25 Fig. 17 shows a final shell element model for 

which the universal file is loaded, and neutral plane 
shell elements are generated from the tetrahedral 
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elements on the basis of the above -described neutral 
plane element division method. The time until the 
final neutral plane shell elements are generated from 
the three-dimensional CAD geometric data is about 30 
min in a PC having a CPU speed of 1 GHz. 

When analysis input data is created by adding 
boundary conditions or analysis conditions 
corresponding to the type of analysis to the geometric 
model, injection molding analysis can be executed. 

As described above, according to the above 
embodiment, a shell element analytical model can be 
automatically created on the basis of three-dimensional 
CAD geometric data, unlike the conventional method of 
manually generating a neutral plane shell element or 
shell element. Accordingly, the time of shell element 
analytical model creation from three-dimensional CAD 
geometric data can largely be shortened. 

As many apparently widely different embodiments 
of the present invention can be made without departing 
from the spirit and scope thereof, it is to be 
understood that the invention is not limited to the 
specific embodiments thereof except as defined in the 
appended claims. 
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